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Topological insulators are a novel materials platform with high applications potential in fields ranging from 
spintronics to quantum computation. In the ongoing scientific effort to demonstrate controlled manipulation 
of their electronic structure by external means - i.e the provision of knobs with which to tune properties - 
stoichiometric variation and surface decoration are two effective approaches that have been followed. In angle 
resolved photoelectron spectroscopy (ARPES) experiments, both approaches are seen to lead to electronic band 
structure changes. Most importantly, such approaches result in variations of the energy position of bulk and 
surface-related features and the creation of two-dimensional electron gases. The data presented here demonstrate 
that a third manipulation handle is accessible by utilizing the amount of super-band-gap light a topological 
insulator surface has been exposed to under typical ARPES experimental conditions. Our results show that 
this new, third, knob acts on an equal footing with stoichiometry and surface decoration as a modifier of the 
electronic band shucture, and that it is in continuous and direct competition with the latter. The data clearly point 
towards surface photovoltage and photo-induced desorption as the physical phenomena behind modifications 
of the electronic band structure under exposure to high-flux photons. We show that the interplay of these 
phenomena can minimize and even eliminate the adsorbate-related surface band bending on typical binary, 
ternary and quaternary Bi-based topological insulators. Including the influence of the sample temperature, these 
data set up a detailed framework for the external control of the electronic band shucture in topological insulator 
compounds in an ARPES setting. Eour external knobs are available: bulk stoichiometry, surface decoration, 
temperature and photon exposure. These knobs can be used in conjunction to fine-tune the band energies near 
the surface and consequently influence the topological properties of the relevant electronic states. 


I. INTRODUCTION 

Topological insulators (TIs) are a realization of a fascinat¬ 
ing and novel state of quantum matter iHHii. Their remarkable 
properties result from metallic, spin-polarized surface states 
which possess linear dispersion: the so-called topological sur¬ 
face states (TSS) Ellll. Synergy between a chiral spin texture 
due to spin-momentum locking and time reversal symmetry 
Ellll means that the TSS are robust to non-magnetic impu¬ 
rities and disorder. Consequently, these systems are generat¬ 
ing interest both in the context of spin generator systems in 
emerging spintronic devices i) and as a platform for the cre¬ 
ation of stable building blocks of quantum information in the 
form of Majorana zero modes ll9i rT0ll . In order to realize such 
exotic applications, it would be beneficial to be able to se¬ 
lectively manipulate the electronic properties of the TSS. As 
viable alternatives to electrostatic gating (e.g. lfTTI - [l4l l. exten¬ 
sive research has been carried out into two different means of 
effective manipulation, namely changes in the bulk stoichiom¬ 
etry of the TI material (e.g. GMIIl) and controlled surface 
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decoration of the TI crystal (e.g. IISMl)- 

The requirement for high bulk resistivity in three di¬ 
mensional TIs has led to a transition in terms of materi¬ 
als from binary teradymites (i.e. Bi 2 Se 3 , Bi 2 Te 3 IIU0) to 
ternary systems (i.e. [Bii_xSbx] 2 Te 3 ll22l |23l), and finally, 
to TlBiSe 2 -based Il24l426l and quaternary Bi 2 Te 2 Se-based 
ifTSl [Thl ITSl l27l compounds (i.e. Tli_xBii+xSe 2 _y ll28l and 
Bi2_xSbxTe3_ySey EtI ISH [30l, respectively). In these sys¬ 
tems, X and y serve as the parameters to fine tune the en¬ 
ergy position of the TSS with respect to the bulk electronic 
states and the Fermi energy. On the other hand, the surface 
engineering route to effective tuning of the electronic band 
structure has been achieved by decorating the surface of a TI 
crystal or film with metals (via deposition) llT9l l20l [3TI43TI . 
via controlled gas adsorption 1201 llU |3lH36l, or by means 
of exposure to the residual gases in a UHV environment (i.e. 
time-dependence in vacuum) EOl [311 [Ml iJTll . 

To enable TIs to reach their full potential, the complete 
parameter space for band structure control and manipula¬ 
tion should be investigated and understood. By contrast to 
the plethora of studies following the two scenarios sketched 
above, only very few reports have been made regarding vari¬ 
ations of the electronic band structure based on other external 
parameters such as temperature variations or exposure to con- 
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tinuous illumination ll^ f38l . 

This paper reports a detailed study of the interplay be¬ 
tween three in-situ factors affecting the surface band structure, 
namely illumination, surface decoration and temperature. The 
analysis considers variations of the bulk conductivity over or¬ 
ders of magnitude via comparison of binary, ternary and qua¬ 
ternary 3D TIs. 

Using Angle-Resolved PhotoElectron Spectroscopy 
(ARPES) at different experimental setups and under variable 
experimental conditions (see Methods), we demonstrate that 
the amount of illumination the surface region has received is 
not a mere secondary factor in determining the surface band 
structure of Bi-based TIs. On the contrary, it will be evident 
that exposure to photons can completely cancel the desired 
modifications induced by deliberate stoichiometric tuning 
of the bulk and/or induced by controlled surface decoration. 
Besides describing the observed photo-induced effects and 
discussing the underlying mechanisms, the investigations 
reported here provide the necessary context in which past 
results on TIs acquired by techniques involving photon 
exposure -such as ARPES- need to be viewed. In doing so, 
these results also provide a framework for future research. 

In the following, we will first present the effects of residual 
gas adsorption and exposure to illumination with a variable 
flux for the case of Bi 2 Se 3 . We then present analogous UHV 
adsorption and photon flux experiments on two, related, qua¬ 
ternary TIs with high bulk resistivity: Bi 1 . 5 Sbo. 5 Te 1 . 7 Se 1 . 3 , 
Bi 1 . 46 Sbo. 54 Te 1 . 7 Se 1 . 3 . The third data section deals with the 
stoichiometric compounds Bi 2 Te 2 Se and BiSbTeSe 2 . The 
main body of the paper finishes with a discussion of the 
underlying physical phenomena and a summary of the main 
consequences of our work for past and future band structure 
studies of these materials. The paper additionally contains 
two appendices in which the effect of sample temperature 
is discussed, for both pristine cleavage surfaces exposed to 
UHV and for intentionally metal-decorated surfaces. 

II. RESULTS: The effect of illumination 
Bi2Se3 

Eig. 1 illustrates our findings on Bi 2 Se 3 , the most widely 
studied topological insulator. In the top panel, we present the 
energy position of (i) the spin degeneracy point of the TSS 
(also known as the Dirac point) and (ii) the onset of the va¬ 
lence band leading edge at fc = — O.llA , both plotted as a 
function of time. In comparison to past studies l|20l|3ll[37l, 
we investigate the effect of an additional parameter which 
turns out to be of great impact: the incident photon flux. Eor 
times up to the point represented in Eig. la by the solid ver¬ 
tical line at ca. 500 minutes, the photon flux was kept at a 
‘low’ value: Fiow= 1.3x10^® photons/(s m^) (see Methods). 
In this regime, both the Dirac point and the valence bands can 
be seen to shift by the same amount to higher binding ener¬ 
gies as time progresses, signalling a rigid shift of the elec¬ 
tronic band structure. This phenomenon has been widely re¬ 
ported as due to downward band bending (BB), as a result 
of a near-surface electrostatic potential UlMlEIllSl. The 
width of the space charge region that accompanies the down¬ 
ward BB has been found to be as large as 20-30 nm from 


ARPES lEOllITllMl and HAXPES studies IM). As the data 
in the right-hand panel of Pig. la show, the time evolution 
of an increasing downward BB can be reversed when - all 
other parameters being held constant - the photon flux in in¬ 
creased by a factor of 100. Upon such strong illumination 
[flux Fhigh=100xE]ow=1.3x 10 ^^ photons/(s m^)], both the 
TSS and near-surface bulk-derived bands shift back to lower 
binding energies, and do so initially at a significantly faster 
rate than the initial, downward shift. Prom the above, it is 
clear that the complete picture in an ARPES experiment is 
highly-dependent on the total photon fluence to which the 
sample is exposed; the result is a dynamic competition be¬ 
tween the downward band bending in UHV and a photo- 
induced flattening of the bent bands. 

According to the most widely accepted scenario, the built- 
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FIG. 1: Time-dependent changes in the electronic hand structure 
of Bi 2 Se 3 under weak and strong illumination (a) Time-dependent 
energy position of the Dirac point (left axis) and the leading edge 
of the VB (right axis). The photon flux is increased by a factor of 
100 at the time marked by the vertical line. Time intervals refer to 
the time since cleavage of the crystal in UHV. The error bars in the 
value of binding energy are estimated to be ±5 meV for both the 
Dirac point and VB. The green solid lines represent single (left panel) 
and double (right panel) exponential functions used to determine the 
time constants (see Eqn. 1) for the downward band bending and the 
fast(slow) upward reversal thereof under high-flux illumination, (b)- 
(d) The near-Ep electronic band structure of Bi 2 Se 3 at different time 
intervals and beam exposures: (b) a freshly cleaved sample without 
any prior exposure to the photon beam [i.e. the first point on curve 
(a)], (c) \(E, k) image acquired after 11 hours from an area of the 
same sample that received only a limited amount of illumination, (d) 
same location on the sample as (b) but now after exposure to both the 
low and high-flux beams [i.e. last point on the curve in the right-hand 
panel of (a)], (e) same location on the sample as (c) after additional 
exposure to higher energy photons (hv = 130 eV) until saturation is 
observed. All data have been acquired at 16K. 
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in electrostatic potential at the origin of the downward BB 
in binary TIs is driven by the decoration of the sample sur¬ 
face with residual gas atoms which act as electron donors 
on chemisorption EOl ISTI l34l [JTll . Strongly-developed BB 
gives rise to a potential well in which also the bulk conduc¬ 
tion and valence band (CB and VB) become quantized lfT9l - 
Izii. Quantum well states in the valence band have been alter¬ 
natively suggested to arise from the expansion of the van der 
Waals spacing between consecutive quintuple layers of the TI 
through the intercalation of gases ll40H42l . In either case, the 
effect of residual gas atoms is central to time-dependent band 
structure modifications in UHV, which are sometimes referred 
to as aging in the ARPES literature. 

An alternative scenario for the downward BB in Bi 2 Se 3 in¬ 
vokes the formation of Se vacancies leading to the formation 
of Se atoms in the gas phase and two extra electrons per va¬ 
cancy in the crystal ll43l . In such a case, one would expect 
that the bands of Se-terminated Bi 2 Se 3 would shift faster than 
Bi 2 Te 2 Se, which is predicted to be Te-terminated ll43l . We 
did not observe a difference in behaviour between these two 
systems in our experiments. Moreover, as discussed later, the 
downward shift is more pronounced at low sample tempera¬ 
tures, which would be somewhat counterintuitive in the Se- 
vacancy scenario. For the above reasons, we also attribute the 
downward BB to the effect of residual gas atoms. 

We observe no or at most marginal differences between the 
time-dependent energy shifts of the surface- and bulk-related 
states (Fig. la). This is a natural consequence of the high 
surface-sensitivity of our ARPES data. Using hi/ = 23 eV, 
we probe the part of the CB and VB localized within 1 nm 
of the surface termination HI; a value which is smaller than 
the out-of-plane spatial confinement of the bulk bands in the 
band-bending potential, and smaller than the spatial extension 
of the surface state in the z-direction, which is 2-3 nm Il45l . 

We attribute the dramatic changes on the sign and rate of the 
observed energy shifts on high-flux illumination to the com¬ 
bined effect of two physical phenomena: 

• surface photovoltage (SPV) effects 

• photo-induced desorption of adsorbed species 

SPV effects are commonly encountered in semiconductor ma¬ 
terials upon illumination with photons of energy exceeding 
the bulk band gap. Light absorption results in the forma¬ 
tion of electron-hole pairs via interband transitions. These are 
separated in the space charge region via the built-in electric 
field that is responsible for the downward band bending. The 
spatial distribution of the SPV-induced charge carriers in the 
band-bending potential opposes the original charge distribu¬ 
tion that gave rise to the band bending. In this way, these 
photo-induced charges cancel the original electrostatic poten¬ 
tial and thus lead to a decrease in the band bending 114614511 . 
Insight into the delicate synergy of SPV effects and photo- 
stimulated desorption of foreign species from the sample sur¬ 
face will be given in Section III. Here we note that surface- 
and bulk-related features shift rigidly to lower binding en¬ 
ergy (i?b) and that the size of the photo-induced energy shift 
is a function of the total photon fluence. The bottom line is 


thus a continuous competition between the adsorbate-related 
charge redistribution and the photo-induced effects; i.e. be¬ 
tween the effects of the residual vacuum (downward BB) and 
of the high-flux illumination (upward shift in energy, neutral¬ 
izing the same downward BB). 

The realization that there is a continuous competition be¬ 
tween adsorbate-related and photo-induced changes clearly 
has consequences for how ARPES can be used to extract 
information on the energy position of the surface and bulk- 
related features in Bi-based 3D TIs. Any such determination 
is incomplete without explicit reference to 

1 . the residual gas pressure 

2 . the time-window the cleaved surface has been in that 
specific vacuum 

3. the photon fluence (i.e. total number of photons re¬ 
ceived per unit area). 

In particular, the last of these three factors has barely been 
addressed in the Bi 2 Se 3 ARPES literature: illumination ef¬ 
fects have been considered in a single study that focused on 
/ii^-dependent changes ll3^ . 

The energy (E) vi. wavevector (k) dispersion for a fresh 
surface of Bi 2 Se 3 (recorded within 35 min. after cleavage) is 
shown in the l(E, k) image of Fig. lb. The Dirac point lies at a 
binding energy of 240 meV - the first red data point in Fig. 1 a 
- and the bulk conduction band rates minimum is located just 
below Ep, although not readily visible for the greyscale used 
(see also Fig. 2d). For such Tow’ flux as the one used for the 
left-hand panel of Fig. la, the effect of residual gas adsorp¬ 
tion prevails over the potential flattening and/or desorption in¬ 
duced by illumination and thus the binding energy of the Dirac 
point increases to 320 meV after 400 min. As can be seen in 
Fig. Ic, at sample positions of the same cleave which have 
not been continuously exposed to the photon beam, the bind¬ 
ing energy of the Dirac point has even higher values. At such 
locations, the BB is now strong enough to act as a potential 
well and clear signatures of confinement are then observed 
in the conduction and valence bands HU |20l |32l. The ob¬ 
servation that unexposed areas are different reveals a degree 
of local character in the photon-induced changes on Bi 2 Se 3 , 
which is dramatically enhanced for bulk-insulating quaternary 
TI compounds 15^ . 

The data shown in Fig. Id correspond to the last data point 
of Fig. la (/ = 630 min). Under higher flux illumination, the 
Dirac point has now shifted back to a binding energy equiv¬ 
alent to that in panel (b). Occupation of the CB states with 
electrons is signalled by the spectral weight close to the Fermi 
level visible inside the linearly dispersing TSS. As the 1{E, k) 
image of Fig. le shows, the band structure can be shifted to 
even lower binding energies via additional exposure to higher 
energy photons (/iz^=130 instead of 23 eV). Comparing pan¬ 
els (d) and (e) with panel (b), one may observe that the ap¬ 
parent energy difference between the CB minimum and the 
Dirac point decreases with increasing hv exposure, while the 
bands themselves appear broadened. These two observations 
are intimately connected as the measured image is a superpo- 
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sition of different areas under the beam-spot, each of which is 
up-shifted by a differing amount. 

The ‘high’ photon flux used to record the (up-shifting) 
data shown in the right-hand panel of Fig. la is of the same 
order as the standard flux that a sample receives at modern 
synchrotron facilities operating in multi-bunch mode with a 
typical storage ring current of 400 mA. The data presented 
in Fig. la were acquired under such conditions at the Swiss 
Light Source, and the same effects have been reproduced at 
the BESSY II source at the Helmholtz Zentrum Berlin. For us 
to observe the well-studied, downward BB and the evolution 
of quantum well states in Bi 2 Se 3 lfT9ll20ll3nf34llJ7]| . the pho¬ 
ton flux had to be deliberately decreased by a factor of 100 , as 
was done for the data in the left-hand panel of Fig. 1 a. With 
the benefit of hindsight, therefore, it is no surprise that the 
pioneering studies on the time-dependent, downward BB of 
TIs were performed with low flux sources, namely laboratory 
setups |l3T][34l [photon flux ^10^® photons/(s m^)] ||53][54l, 
low-flux synchrotron sources lElEl or large-scale facilities 
operated at low storage ring current such as is the case in 
so-called single-bunch mode ll20l . Under standard, high flux 
conditions (400 mA ring current; insertion device beamline), 
the illumination effects are so strong at low temperature that 
they can essentially prevent the time-dependent, downward 
BB. 


Quaternary compounds with high bulk resistivity 

In bulk crystals, the transport properties of the topological 
surface states in Bi 2 Se 3 and Bi 2 Te 3 are difficult to investi¬ 
gate, as bulk conduction due to intrinsic impurities and crys¬ 
tallographic defects dominates. Although progress has been 
achieved on charge carrier doping Il55ll56ll . thin film engineer¬ 
ing and electrostatic gating CHEl, an alternative has been 
explored in the development of bulk materials with reduced 
bulk conduction. The quaternary system Bi 2 _xSbxTe 3 _ySey 
(BSTS) compound HT) has been at the forefront of these 
efforts. We have also investigated the effect of low- and 
high-flux illumination on crystals of two BSTS compositions: 
Bi 1 . 5 Sbo. 5 Te 1 . 7 Se 1,3 (BSTS 1.5) and (BSTS 1.46). In contrast 
to Bi 2 Se 3 , both BSTS1.5 and BSTS1.46 exhibit insulator-like 
resistivity curves with typical low temperature values of or¬ 
der 10 Ohm cm (Fig. 2a). The low-tempearature resistivity of 
these quaternary compounds is more than 4 orders of magni¬ 
tude higher than the resistivity of Bi 2 Se 3 . 

Low flux regime 

The binding energy of the Dirac point for freshly cleaved 
BSTS 1.46 and BSTS 1.5 samples typically ranges between 0.1 
and 0.2 eV (Fig. 2). This uncertainty arises from the fact 
that the downward BB sets in at its maximum rate in the first 
minutes after cleavage (Fig. 2c). Unlike the case of Bi 2 Se 3 , 
cleaved BSTS 1.5 and BSTS 1.46 samples are characterized by 
the complete absence of a bulk conduction band at the Fermi 
level, as seen in the near-F^p energy dispersion (Fig. 2b, after 
0.6 h in UHV). The comparison shown in Fig. 2d between the 
fc|| = 0 energy distribution curves for freshly cleaved Bi 2 Se 3 
and BSTS 1.5 (the former data trace is a vertical slice through 
the data of Fig. lb), showing that electron occupy the conduc- 
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FIG. 2: The effect of low-flux illumination on quaternary TIs 
with high bnik resistivity: Bii.sSbo.sTei.TSei.a (BSTS1.5) and 

Bii, 46 Sbo. 54 Tei, 7 Sei ,3 (BSTS1.46) (a) T-dependent resistivity of 
BSTS1.5, BSTS1.46 and Bi 2 Se 3 . (b) The near-Ep electronic band 
structure of BSTS 1.46 after only minimal exposure to photons and 
UHV exposure as indicated, (c) BSTS 1.5: Time-dependent energy 
position of the Dirac point (left axis) and the intensity maximum of 
the fc-integrated valence band spectrum within the energy window 
of measurement (right axis) under weak illumination [1.3x 10^® /(s 
m^)]. Error bars in the value of binding energy are estimated to be 
±10 meV for the Dirac point and ±5 meV for the valence band. The 
green solid line represents an exponential fit to the time-dependence 
of the downward band bending of the Dirac point (Eqn.l). (d) Com¬ 
parison of the energy distribution curves at fc|| = 0 for freshly 
cleaved BSTS 1.5 and Bi 2 Se 3 (both 35-40 min. after cleavage), re¬ 
vealing the absence of a low binding energy peak signalling the con¬ 
duction band in BSTS. The Bi 2 Se 3 spectrum is from the data shown 
in Fig. lb. All ARPES data were acquired at 16K from single crystals 
from the same batch as those used for the resistivity measurements. 


tion band states only in the binary TI. 

However, downward band bending in BSTS can also bring 
the bulk conduction band well below Ep, as is seen in the 
right-hand panel of Fig. 2b for BSTS 1.46. By tracking the 
complete time dependence (Fig. 2c) of the surface band struc¬ 
ture of BSTS 1.5 at low photon flux (flux Flow), the first signs 
of CB occupation are seen within a time-scale of 100 min for 
these vacuum conditions (i.e. P = 5 x 10“^^ mbar) at low 
temperature (T = 16 K). The data shown in Figs. 2b and 2c 
are therefore a vivid demonstration that also in BSTS-based 
compounds, time-dependent BB should also be taken into ac¬ 
count when discussing the energy position of the Dirac point 
with respect to Ep. The data shown in Eigs. 1 and 2 illus¬ 
trate that variation in the time interval between cleavage and 
measurement can yield changes in the observed binding en¬ 
ergies that are greater than the intentional changes made via 
alteration of the chemical composition of the bulk crystals. 

A comparison of the early data points in Eigs. la and 2b re- 
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veals that the downward band bending evolves faster in BSTS 
than in Bi 2 Se 3 under equivalent conditions. This means that 
while BSTS reaches a saturated energy shift after approxi¬ 
mately 150 mins of UHV exposure to a low photon flux, the 
analogous, downward energy shift in Bi 2 Se 3 fails to saturate 
within the time window of our measurement. The saturation 
energy value for the BB in BSTS shown in Fig. 2c represents 
the outcome of a continuous competition between decoration 
with residual gas atoms and exposure to the low-flux photon 
beam. Just as in the case of Bi 2 Se 3 (Fig. Ic), at sample lo¬ 
cations that have been only minimally exposed to photons, 
greater band bending is observed. 

In order to provide a quantification of the energy shifts, we 
have fitted an exponential function to the time-dependent data 
shown in Figs, la and 2c: 

E{t) = ±Eoe-^/^ (1) 

with a positive prefactor describing the downward BB (in 
which the Dirac point binding energy increases with time) 
and a negative prefactor, the upward shift under high flux 
conditions. The characteristic time constant r for the band 
bending is 193 min for Bi 2 Se 3 and 32 min for BSTS, the latter 
echoing the ‘steeper’ temporal evolution of the downward 
shift in BSTS. A possible origin for this behaviour could 
lie in the more effective screening of the adsorbate-induced 
surface charges underlying the downward BB in Bi 2 Se 3 , 
compared to the bulk-insulating BSTS. We also note that 
as the chemical composition of the termination layer of the 
two compound differs, this could also impact the sticking 
coefficient for residual gas atoms, and therefore also modify 
the BB behavior. 


High flux regime 

Fig. 3 shows the influence of a high-flux photon beam on 
BSTS. High flux means Fhigh= 1.3x10^^ photons/(s m^). 
The starting point is a saturated downward BB surface. Un¬ 
der this strong illumination, the binding energy of the Dirac 
point is seen to reduce, indicating that the bands flatten and 
the overall band structure is shifting upwards to lower binding 
energy (Fig. 3a). As was observed for Bi 2 Se 3 (right panel of 
Fig. la), this photon-induced effect is rather rapid in the first 
few minutes of high flux exposure. 

The 1{E, k) images in panels (b), (c) and (d) of Fig. 3 
show the observed band structure changes in BSTS on in¬ 
creasing high-flux beam exposure, in which pronounced spec¬ 
tral broadening is apparent, taking the form of a shadow¬ 
like feature sutTounding the main experimental bands. As in 
Bi 2 Se 3 , the broadening is a manifestation of the local action 
of the photon beam in which distinct electronic band struc¬ 
tures are summed, each up-shifted differently in energy due 
to the different number of photons in the spatially inhomoge¬ 
neous profile of the photon beam. The data of Fig. 3 show 
the local character of these photon-induced effects to be fur¬ 
ther enhanced in the case of bulk-insulating BSTS. This effect 
has been modelled in Ref. |52l and exploited to create micro¬ 
metric patterns in the topological band structure. 

Unlike the downward BB in Figs, la (left panel) and 2c, the 
time evolution of the upward energy shift on high-flux illumi- 
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FIG. 3: The effect of high-flux illumination on quaternary TIs 
with high hulk resistivity: Bii.sSbo.sTei.rSei.s (BSTS1.5) and 

Bii, 46 Sho. 54 Tei, 7 Sei ,3 (BSTS1.46) (a) BSTS1.46: Time-dependent 
energy position of the Dirac point (open circles) under the influence 
of strong illumination (yellow panel) and after continuous illumi¬ 
nation ceased (grey panel). The energy position of the Dirac point 
recovers only partially once the light is off, and does so on a much 
longer time-scale than the onset of the photo-induced upward shift 
in the yellow panel. Green solid lines are exponential fits on the 
data points using Eqn. 1. As described in the text, the ‘illuminated’ 
region is best fitted with a double exponential function whose con¬ 
stituent curves are shown by pink and blue dashed lines. A single 
exponential fit can accurately describe the ‘dark’ region. Arrows in 
panel (a) [i.e. BSTS 1.46] denote the time intervals where the anal¬ 
ogous data on BSTS 1.5 [panels (b) - (d)] has been acquired. The 
error bars in the value of binding energy are between ±5 and ±15 
meV. (b) BSTS 1.5: near-Ep electronic dispersion upon minimal pho¬ 
ton exposure after considerable UHV exposure (i.e. 10 hours after 
cleavage), (c) Same as (b) but after 15 min high-flux photon expo¬ 
sure [1.3x 10^^ photons/(s m^)]. (d) Same as (c) but with the total 
high-flux exposure time now increased to 2 hours. All data have been 
acquired at 16K. 


nation cannot be accurately fitted with a single exponential. 
For both Bi 2 Se 3 and BSTS a double exponential function is 
required. The resulting fits are shown using solid green lines 
in the left-hand panel of Fig. 3a for BSTS and in the right- 
hand panel of Fig. la for Bi 2 Se 3 . For both materials, the dou¬ 
ble exponential fits contain a ‘fast’ time constant (rp) that de¬ 
scribes the first few minutes of the upward shift, and a ‘slow’ 
time constant (T 2 ) that tracks the evolution at later stages. For 
the BSTS data shown in Fig. 3, the ti(t 2 ) curves are shown in 
Fig. 3a using pink(blue) dashed lines. In Section III, we bring 
together all the exponent fit-data in Table 1, and discuss how 
the r values in the double-exponents are consistent with the 
combined effects of surface photovoltage (faster) and slower 
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photo-induced desorption processes. 

The final step in the presentation of the ARPES data from 
BSTS is the temporal decay of the upward SPV-induced en¬ 
ergy shift - i.e. the return to a time-depedent increase in the 
binding energies of the Dirac point and other characteristic 
features of the electronic structure - when the high-flux illu¬ 
mination source is turned off. The right panel of Fig. 3a shows 
that on high-flux switch-off, the downward BB begins to re¬ 
assert itself, and the topological surface state returns to higher 
binding energy, a process which can be captured by a third 
exponential relation shown in the figure by the solid, green 
line. 

In analogy to the situation in conventional semiconductors, 
we attribute the sluggish nature of this return to downward 
BB to the very small drift velocity across the space charge 
region (SCR), meaning that photo-induced charges can exist 
well beyond the exposure time, depending on temperature 
and on the magnitude of the original band bending. Lifetimes 
of photo-induced charges of up to 10^ seconds have been 
reported for GaAs ll57l . 

Other ternary and quaternary compounds 

Before Section III presents a discussion of the relative mer¬ 
its of SPV and photon-stimulated desorption as mechanisms 



FIG. 4: The effect of high-flux illumination on Bi 2 Te 2 Se (BTS) 
and BiShTeSe 2 (BSTS2) (a) T-dependence of the resistivity of 
BSTS2 (both ‘type P and ‘type 2’ [for definitions, see text]) and 
BTS. (b) Time-dependent energy position of the Dirac point of type 
1 BSTS2 under the influence of a high-flux photon beam, with data 
from the valence band edge superimposed as described in the main 
text (c) BSTS2: near-Sp electronic band structure for a type 1 sample 
at times as indicated by the arrows in panel (b). (d) Time-dependent 
energy position of the Dirac points of BSTS2 (‘type 2’, blue sym¬ 
bols) and BTS (black symbols) under the influence of a high-flux 
photon beam. The grey area denotes a time-period in the BTS ex¬ 
periment that was without high-flux photon exposure. All data have 
been acquired at 16K. Resistivity and ARPES measurements were 
performed on the same crystal pieces. 


behind the reversal of the downward BB under high-flux il¬ 
lumination, the data collection feeding this discussion will be 
completed with data from other ternary and quaternary 3D TI 
compounds. 

Fig. 4 summarizes the data for the two nominally stoichio¬ 
metric compounds Bi 2 Te 2 Se (BTS) and BiSbTeSe 2 (BSTS2). 
Panel (a) shows the temperature-dependence of the resistivity 
for these compounds. BTS shows a similar resistivity curve to 
the two BSTS compositions studied in Figs. 2 and 3, and the 
low-T saturation value of the resistivity is only 20% smaller 
than for BSTS 1.46. 

Resistivity curves of BSTS2 were found to vary signif¬ 
icantly from sample to sample. More than 20 samples 
have been measured and two characteristic cases which were 
widely encountered are shown in Fig. 4a. BSTS2 samples of 
what we refer to here as ‘type 1 ’ (red trace in Fig. 4a) retain 
an insulating up-turn in resistivity, but the low-T value is 2 
orders of magnitude lower than for the BSTS samples (Fig. 
2a) and BTS (Fig. 4a, black trace). Moreover, for ‘type 1’ 
BSTS2 samples the resistivity exhibits a clear transition to a 
more metallic behavior at higher temperatures. The overall, 
temperature-dependent resistivity curve of ‘type 1’ BSTS2 is 
reminiscent of that in extrinsic semiconductors in which, as 
the temperature is raised, the resistivity initially decreases due 
to the thermal excitation of electrons trapped at impurity sites 
followed by a high-temperature increase as the carrier mobil¬ 
ity decreases due to thermal fluctuations. In short, we attribute 
differences in the resistivity curve of BSTS2 samples of ‘type 
1 ’ and other ternary and quaternary compounds to the larger 
bulk defect density of the former. On the other hand, BSTS2 
samples we refer to as ‘type 2’ (blue trace in Fig. 4a) show 
an insulator-like T-dependence of the resistivity at all temper¬ 
atures. The low-T resistivity value of ‘type 2’ BSTS2 sam¬ 
ples is very similar to BTS (Fig. 4a) and other BSTS samples 
(Fig. 2a). The ARPES signatures of the two different types of 
BSTS2 samples after exposure to a high-flux photon beam are 
very instructive for establishing a link between the bulk car¬ 
rier concentration and their response to external illumination 
in this system. 

Figs. 4b-4d present our ARPES results on BTS and BSTS2 
after exposure to illumination comparable to the high-flux 
beam used in Figures 1 and 3. The electronic structure of 
‘type 1’ BSTS2 samples is very sensitive to external illumina¬ 
tion. Fig. 4b shows that the Dirac point was shifted upward 
to the Fermi level after only 17 minutes of high-flux expo¬ 
sure. To be able to determine the trend at later times, as a first 
step, an estimate of the density of states is made by integrat¬ 
ing the 1{E, k) images along the momentum direction. The 
shift of the valence band edge in these effective DOS data is 
then tracked as a function of time and the resulting binding 
energy vs. time curve is re-scaled, so it matches that from the 
Dirac point over their common range. The data of panel (b) 
show a turnover to a slower evolution shortly after the Dirac 
point crosses Ep- Although the data do not show a satura¬ 
tion of the upward shift, it is clear from the I(i?, k) images of 
‘type 1’ BSTS2 shown in Fig. 4c that for the highest expo¬ 
sures measured, the Fermi level is definitely below the Dirac 
point, pointing to the p-type character of this material. Fig. 
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4d shows that the BSTS2 samples with large low-T resistivity 
(‘type 2’) also show an upshift of the Dirac point upon high 
flux illumination, but this takes place slower than for the ‘type 
r crystals, and in the time window for which there are data, 
there are no signs of a slowing or saturation of the up-shift. 

The overall behavior of BTS on increasing photon fluence 
(Fig. 4d black markers) is similar to ‘type 1’ BSTS2 shown in 
panel (b), but the magnitude and rate of the upward shift are 
much smaller in BTS. The shaded area in panel (d) denotes 
a time period for which the photon beam illuminated other 
(neighboring) locations on the BTS sample. During this time 
we estimate that only a small amount of radiation affected the 
sample location from which these ARPES data points have 
been obtained. Interestingly, the rate of the initial upward 
shift in BTS - i.e. that which would be governed by a time- 
constant Ti - is very similar to that of ‘type 2’ BSTS2 (Fig. 
4d empty blue symbols). As none of the binding energy vs/ 
time curves in Fig. 4 saturate, we refrain from carrying out 
(double-)exponential fits of the temporal behavior. 

Summarizing the main points arising from the data pre¬ 
sented so far, it is clear that competing trends in the energy 
position of the topological surface states, as well as the 
valence and conduction bands in Bi-based 3D topological 
insulators result from UHV exposure (linked to adventitious 
adatom adsorption) and high-flux photon beam exposure. 
This implies that the energy position of the near-surface 
electronic band structure of different TI compounds can be 
fine-tuned, after the bulk stoichiometry of the bulk crystal has 
been determined in the growth process. In the next section, a 
discussion of the possible mechanisms behind the observed 
changes is presented, followed by examination of insight that 
can be gained by considering the sample temperature as a 
further parameter. 

III. DISCUSSION 

Two control ‘knobs’ 

In the above, we have identified and characterized two de¬ 
terministic, controllable and (given long timescales) reversible 
parameters (or ‘knobs’), each of them operative after the syn¬ 
thesis of the topological insulator material in question. The 
first control parameter is the well-known, adsorbate-induced 
band bending 11211 [34l |431 . which leads to an increase of the 
Dirac point binding energy. An alternative interpretation of 
these observations could be band bending due to the intrinsic 
properties of the surface of what is an extrinsic (doped) semi¬ 
conductor. However, for n-type Bi 2 Se 3 Il43ll58ll . this mecha¬ 
nism would predict upward BB Il46ll47l . the reverse of what 
is observed experimentally ifT^llOl . Therefore the assignment 
to adsorbate-induced BB is the more convincing of the two. 

A second external control parameter is based on the inter¬ 
action of the near-surface region of the topological insulators 
with a photon beam. This interaction leads to a decrease of 
the Dirac point binding energy, the energy value of which can 
be tuned via choice of the photon flux, energy and exposure 
time. This second knob operates at a local spatial level ll52l . 


The physical phenomena induced by high-flux illumination 

The data suggest that the underlying mechanism behind the 
response to high-flux illumination is a combination of the sur¬ 
face photovoltage (SPV) effect and of the photo-induced des¬ 
orption of adsorbed atoms. Before discussing the delicate in¬ 
terplay between these factors in more detail, we discuss some 
alternative mechanisms. 

• Negligible local heating from the high-flux beam. Cal¬ 
culation of the local, beam-induced temperature rise of 
the sample for a 23 eV photon beam with an output of 
1.3x10^^ photons/(s m^) and with a Gaussian profile 
of FWHM = 30 /rm (typical experimental conditions), 
arrives at maximally a few mK 1591 . This effectively 
excludes beam-induced heating as the origin of the ob¬ 
served upward shifts of the bands. 

• Photo-induced creation of vacancies is not the dom¬ 
inant process. We observe the same kind of photo- 
induced, upward shifts on Se-terminated Bi 2 Se 3 , as for 
Te-terminated Bi 2 Te 2 Se, and for systems with mixed 
Se-Te termination such as BSTS l43l . This observation 
cannot easily be consistent with the creation of Se (or 
Te) vacancies on the surface of TIs after exposing them 
to high-flux photon beams. 

Apart from the above mentioned mismatches for alterna¬ 
tive scenarios, the observation of two different time-scales for 
the upwards energy shifts (fitted using double exponentials) 
would seem to fit to an interplay of SPV and photo-induced 
desorption phenomena. SPV is an electronic effect and can 
be reasonably expected to be faster than the desorption of 
adatoms. Since the ti values extracted (see Tab. 1) describe 
processes some twenty times faster than T 2 , it would be natu¬ 
ral to associate ti with SPV and T 2 with photo-induced des¬ 
orption. However, it should be stressed that a combination of 
SPV and desorption is necessary to describe the reported data 
of Figs. 1-4. In particular: 

• The ‘few-hour’ timescale observed in Figs. 3 and 4 
for the saturation of the photo-induced upward shifts 
goes well beyond anything reported for surface photo¬ 
voltage effects in regular semiconductors such as Si 
or GaAs. In these systems, SPV has been shown to 
saturate on time-scales ranging between only picosec¬ 
onds to microseconds l57ll60l . so the long-timescle ef¬ 
fects seen for these TI surfaces can also involve photon- 
stimulated desorption. Having said that, we do note that 
SPV time-scales of a few minutes, comparable to the 
Ti’s found in this study, have been reported for the SPV 
onset in ZnO nanorods M- 

• The sensitivity of the photo-induced upward shift to 
the photon energy (see Fig. Id vi. le: greater for 
higher photon energy) does argue for a role for the SPV- 
mechanism and is in keeping with a recent report on 
Bi2Se3 ESI. 

• The temperature dependence of the energy shifts pre¬ 
sented in Appendix 1 of this paper show qualitative 
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agreement with the expectations for a process in which 
SPV plays a significant role. 

One aspect in which the observed behavior differs from 
the standard SPV expectation regards the Fermi level. 
SPV-induced charges accumulating at the surface of a semi¬ 
conductor should lead to the deviation of the Fermi level from 
the equilibrium position ll46l . yet we do not observe this. This 
suggests that also for ternary and quaternary Bi-bassed 3D 
TIs the topological surface states are important in maintaining 
charge neutrality at the surface, as they have been argued to 
be in Bi 2 Se 3 ll38l . A difference between the bulk-insulating 
systems and Bi 2 Se 3 (with its degenerate, defect-doped bulk) 
is that it is possible that photo-induced charges can be trapped 
locally in the sub-surface region. When swept apart by the 
in-built potential which gives rise to the downward band 
bending, these charges may therefore alter the energetics 
of the surface electronic structure in a manner somewhat 
analogous to a local back-gate, without the need for explicit 
surface charging and the requirement of a measurable shift of 
the Fermi level at the surface from the equilibrium position 

IMl. 

The challenge of ohtaining flat hand conditions 

The data of Figs. 1 and 2 show that low-flux ARPES spectra 
acquired on binary, ternary and quaternary Bi-based TIs cor¬ 
respond to a band-bent situation E0ll2ni3T1[^[37l . as the ini¬ 
tial, downward shift evolves very quickly in the first few min¬ 
utes after cleavage. It is therefore questionable whether flat 
band conditions can be captured experimentally within this 
short time window after cleavage. Although the data of Figs. 
1-4 suggest that photo-induced effects could be used to neu¬ 
tralize the downward band bending, before any conclusions 
are drawn as regards the true, flat-band electronic structure of 
the TI systems studied, it is appropriate to discuss briefly the 
differences between considering an extrapolation of low-flux 
or high-flux data to the limit of infinite fluence; 

• Low flux case. If the exponential fits made to model the 
time dependence of the Dirac point energy - for exam¬ 
ple in the left-hand panel of Fig. la or in Fig. 2c - were 
used to extrapolate to ‘infinite fluence’, it is clear that 
the Flow case cannot yield flat band conditions. This un¬ 
derlines that there is more afoot than only the creation 
of (very) long-lived electron hole pairs in the reversal of 
the downward band bending. 

• High flux. The time dependences observed in Figs, la 
(right). Fig. 3a (left) and Fig. 4b for Bi 2 Se 3 , BSTS and 
BSTS2 show that the extrapolation to ‘infinite fluence’ 
using Fhigh conditions yields a principally different end¬ 
point than the Flow case. The high-flux limit is clearly 
much closer to flat band conditions than the low-flux 
case. 


Taking this point to heart. Fig. 5 illustrates how the best 
achievable approximation of flat band conditions can be gen¬ 
erated using external, high flux illumination. The vertical 



photon flux X excitation energy [10^^ photons/(m^s) x eV] 


FIG. 5: The effect of photon flux and excitation energy in ap¬ 
proaching flat band conditions The vertical axis is the minimum 
binding energy (Es) of the Dirac point, arrived at by infinite ex¬ 
trapolation of exponential fits to the time dependence on high-flux 
illumination. Symbols represent the Eb values for different 3D TI 
compounds. The left-most data points underlaid in green show the 
maximally band-bent starting situation before high-flux illumination. 
The top axis gives the photon fluxes; and the bottom axis the prod¬ 
uct of photon flux and photon energy. The orange guides to the eye 
show how the Dirac point binding energy decreases as the flux is 
increased and the purple guides to the eye illustrate further band flat¬ 
tening when, already at high flux, the photon energy (and hence the 
product flux X hiy) is increased. The right-most data-point of both 
BSTS samples and the BSTS2 sample lie on the light blue back¬ 
ground, and thus, from a close-to-flat-band ARPES perspective can 
be considered to be p-type. On the other hand, for Bi 2 Se 3 and BTS 
the analogous Dirac point Eb remains on the grey background, mak¬ 
ing these materials n-type. 


axis shows the binding energy of the Dirac point obtained 
by extrapolation of the double exponential functions that de¬ 
scribe the upward shift of the bands (see for example Fig. 3a 
[yellow background]) to infinite time. The starting point is the 
left-most group of symbols (superimposed on a green back¬ 
ground). These are not extrapolated values, but are the Dirac 
point energies from I(F, k) images showing the most down¬ 
wardly band bent situations encountered, prior to high-flux 
illumination. The first step towards flat-band conditions is 
through increasing the photon flux as shown in the top hor¬ 
izontal scale. This leads to an evolution of the extrapolated 
Dirac point binding energies as shown by the orange guides,. 

For the BSTS2 sample (solid black diamond in Fig. 5), this 
is already enough to put the extrapolated Dirac point more 
than 100 meV above the Fermi level, clearly indicating this 
material to be - in the absence of downward band-bending - 
p-type. For the other materials, with maintenance of the high 
flux, the photon energy was increased. The bottom axis shows 
the product of the flux and the photon energy, and this num¬ 
ber can be seen to represent the ability of the high-flux photon 
beam to create (multiple) electron-hole pairs, freeing charges 
to flatten the band bending at the surface of the crystal. The 
extrapolated Dirac point energies of the BSTS samples (green 
triangle, red filled circle) then shift further to still lower bind¬ 
ing energy (see purple guides), finally landing just above the 
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TABLE I: From band bent (BB) to flat band conditions - compilation of key data for various 3D topological insulators. This table brings 
together ARPES data on the electronic structure of a number of Bi-based TI systems after downward BB, and on how flat-band conditions for 
these systems can be generated using high-flux photon exposure. These data are compared to the low-T electrical resistivity of each system. 
Row 1: the initial binding energy of the Dirac point (J5d) is extracted from the figure panels indicated which track the evolution of the upward 
energy shift on high-flux illumination. Row 2: time constants for the double-exponential fit to the high-flux-induced upward shift of the 
electronic structure to lower binding energy. Row 3: flat-band Eo values. Minimum binding energy of the Dirac point, from extrapolation 
of the exponential fits to the time dependence when the product flux x/it/ has its highest value for each compound (i.e. right-most data points 
for each symbol in Fig. 5). Row 4: this 5 minute time window for up-shifting corresponds to that described by ri, where double-exponential 
fits have been made. Shown is the upshift as a percentage of the difference between the values shown in Rows 1 and 3. High degrees of BB- 
reversal are encountered at short exposure times in samples showing relatively low resistivity. Row 5 p at T=16K. Row 6: for the designation 
n|p-type, the Fermi level - under flat band conditions - is located at the edge of the bulk CB| VB, as derived from the data of Row 3 and the 
I(i?, fc) images from ARPES. Row 7: largest observed Eo values, presumably under saturated (downward) band bent conditions (see also Fig. 
6 ). Row 8: gives the maximum value of downward BB observed on exposure to residual UHV gases, given by the difference between the data 
of Row 7 and the flat-band E^ values of Row 3. 


No. 

Quantity 

Bi2Se3 

BiSbTeSe 2 
(type 1) 

Bii,46Sbo.54Tei.7Sei.3 

Bi2Te2Se 

BiSbTeSe 2 

(type2) 

1 . 

Initial Ed (meV binding energy) 

325 

235 

325 

360 

290 


This paper, figure panel 

la 

4b 

3a 

4d 

4d 

2 . 

Ti, T 2 (both in min.), on exposure to TTigh 

0.6, 14.8 


4.3,81.8 



3. 

Ed (meV binding energy) at flat bands 

160 

-110 

-10 

185 

-110 

4. 

% of initial, downward BB reversed after 5 min. at Thigh 

69% 

52% 

13% 

1.4% 

0 .6% 

5. 

Electrical resistivity (T=16K, in Ohm-cm) 

<0.001 

<0.1 

11 

8 

8 

6 . 

n- or p-type at flat band conditions? 

n-type 

p-type 

p-type 

n-type 

p-type 

7. 

Max. Ed (meV binding energy), via UHV exposure 

440 

290 

360 

470 

290 

8 . 

Max. BB (meV, downward), via UHV exposure 

280 

400 

370 

285 

400 


Fermi level for the highest value of the product of flux and 
photon energy. This suggests that like the case of BSTS2, the 
surface region of the BSTS-based materials would be p-type, 
in the absence of downward band bending. In contrast, the 
extrapolated Dirac point binding energies of Bi 2 Se 3 (blue tri¬ 
angle) and BTS (white circle), remain steadfastly more than 
150 meV below the Fermi level, indicating these systems are 
- also without downward bend bending - n-type. 

The combination of high photon flux and higher photon 
energy that yields the right-most symbols for the Dirac 
point energy for each TI compound in Fig. 5 represents as 
close to flat band conditions as can be practically achieved 
in an ARPES experiment. Table 1 gathers together the 
key electronic structure data from the different materials 
investigated in this manner. Even though the strategy adopted 
seems successful in removing the distorting effects of the 
near-surface downward band-bending, it should still be recog¬ 
nized that the real-life situation is still a dynamic equilibrium 
between the upward-shifting effects of the photon flux (SPV, 
photon-stimulated desorption) and the downward-shifting 
adsorbate-induced band-bending. 

Comparison of the response of the various TI compounds to 
high-flux illumination 

In Section II, it was described how high-flux illumination 
can reverse adsorbate-related downward BB and in Section 
III the pro’s and con’s of the different physical phenomena 
were reviewed in favor of a combination of surface photo volt¬ 


age and desorption effects. Table 1 brings together not only 
the key electronic structure data for the maximally bend bent 
condition but also those relevant to the ‘flat band’ situation. 
In addition. Tab. 1 shows the rate at which the band bent sit¬ 
uation was reversed, expressing the up-shift within the first 
five minutes of high-flux illumination in terms of a percent¬ 
age of the original (downward) band bending. In particular, it 
is clear that the crystals with low resistivity at the cryogenic 
temperatures at which the ARPES data in the main body of the 
paper were recorded - and that means Bi 2 Se 3 and the ‘type 1’ 
crystals of BSTS2 - exhibit a fast upshift of between 50 and 
70 % of the original band bending. In particular, comparison 
of the two types of BSTS2 crystal is interesting, whereby the 
fast upshift of the ‘type 2’ samples (which have a high resis¬ 
tance at low T) is less than 1% of the original downward band 
bending. 

Assuming a simple Schottky model for which the maximal 
BB exceeds the thermal broadening, the solution of Poisson’s 
equation provides a relationship between the band bending 
potential {V) and the width of the space charge region {d): 


eeo 

where A), is the bulk carrier density, e is the static dielectric 
constant and Eq is the vacuum permittivity 1 ^ . 

If the bulk carrier density increases, as is the case on going 
from the ‘type 2’ to the ‘type 1’ BSTS2 crystals, one would 
expect the the width of the space charge region (SCR) to de- 
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FIG. 6: From saturated UHV-induced band bending to flat band conditions by bigb-flux ilinmination. (Top row) Experimental l(E, k) 
images arrived at by exposure of different compounds to residual gases in UHV (P < 2 x 10“^°mbar) for time periods in the range of 10-48 
hours, during which external illumination of these specific sample locations was kept to a minimum and the downward band bending was seen 
to saturate (see also row 7 of Table 1). (Bottom row) Experimental 1(E, k) images recorded after high-flux, high photon-energy exposure, 
ensuring close to flat-band conditions. The energy of the Dirac point for each data panel in the lower row differs by less than 20 meV from 
last (right-most) data point for each corresponding material plotted in Eig. 5 (see also row 3 of Table 1). Samples whose data are on a grey 
(light blue) background correspond to «-type (p-type) compounds and their Fermi energies are pinned within a few meV of the bottom of the 
conduction band (top of the valence band). The four schematic illustrations show the electronic band structure for band bending (top row) and 
flat-band (bottom row) conditions, for both «-type (left) and p-type (right) compounds. Data have been acquired at temperatures between 16K 
and 35K, and all \(E, k) images are plotted on the same axes and scales as Bi 2 Se 3 . 


crease for a given total band bending potential. By means of 
parametric solutions of all equations related to the band bend¬ 
ing and charge transfer, recently Brahlek et al. provided con¬ 
tour plots showing the relation between the SCR and N^uik 
for 3D topological insulators ll63l . 

If the SCR for the ‘type 2’ BSTS2 and the other highly re¬ 
sistive Bi-based TIs were indeed broader, this could slow the 
carrier diffusion across the SCR, and in turn this would slow 
the charge redistribution leading to the flattening of the down¬ 
ward BB - regardless of whether SPV or stimulated desorption 
of adsorbates is at work. 

Within this framework, the width of the space charge re¬ 
gion in BSTS1.46 (Ni,=2 x 10^® cm“^) can be expected to 
be more than 10 times larger than for Bi 2 Se 3 (iVf,=5 x 10^^ 
cm“^), taking the bulk carrier concentrations from Ref. 
Although this model is too crude to capture the details of the 
adsorbate-induced band bending, it provides important insight 
which is also in line with the results presented here. The trans¬ 
port data tell us that Bi 2 Se 3 has the most bulk carriers, fol¬ 
lowed by ‘type 1’ BiSbTeSe 2 . Correspondingly, these are the 
two Bi-based TI compounds we studied that exhibit the fastest 
rates for the reversal of the band bending (see Figs, la and 4b), 
which would be in keeping with their possession of a narrower 
SCR. In contrast, the more resistive Bii. 46 Sbo, 54 Tei. 7 Sei. 3 , 


Bi 1 . 5 Sbo. 5 Te 1 . 7 Se 1 . 3 , ‘type 2’ BiSbTeSe 2 and Bi 2 Te 2 Se have 
a broader space charge region and display more sluggish BB 
reversal rates on high photon flux exposure. 

An independent indication of the relative extension of the 
SCR is given in Fig. Ic and the right panel of Fig. 2b. Clear 
confinement effects are seen in the electronic band structure of 
Bi 2 Se 3 after 11 h UHV exposure (Fig. Ic), while no signs of 
quasi-2D confinement is seen at all in the measured electronic 
structure of BSTS1.46, even after exposure to the same UHV 
conditions for 48 hours (Fig. 2b, right). Again, this speaks for 
a narrower SCR in Bi 2 Se 3 , compared to BSTS. 

The energy position of the Fermi level deep inside the bulk 
of a 3D topological insulator is determined by its bulk carrier 
density, Ni, Il63l . When Nt, is smaller than 10^® cm“®, Ep is 
pinned close to the conduction band minimum for an n-type 
material or similarly close to the valence band maximum for 
a p-type system ll63l . For the surface region of either type 
of material - and it should be stressed that it is this region of 
the sample that is relevant for ARPES, STM/STS or transport 
experiments for true bulk insulating samples - the same holds 
only if flat band conditions can be generated. Consequently, 
given that the high-flux, high photon-energy ARPES data are 
as close as one can come to flat band conditions in the context 
of a surface sensitive experiment (Table 1, row 3), the spec- 
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troscopic data presented here can be used to sort the Bi-based, 
3D TI compounds we have studied into either n-type or p-type 
electronic structures in the near-surface region. This is done 
in the data shown in Table 1. 

To bring home this important point in a graphical manner, 
Figure 6 shows - for each of five different Bi-based 3D topo¬ 
logical insulators studied - the juxtaposition of the l(E, k) im¬ 
age for the maximally downward BB ‘starting’ situation, with 
the analogous l(E, k) image for the most up-shifted situation 
observed, recorded using high-flux illumination. The top row 
of Fig. 6 is also the state of affairs one would expect to be 
confronted with in the moderate vacuum encountered in cryo¬ 
genic apparatus used to measure transport properties of bulk 
crystals. Under such conditions, the conduction bands of all 
of the Bi-based 3D TIs we have investigated are partially oc¬ 
cupied - they possess n-type surfaces, with topologically triv¬ 
ial states mixed in with the TSS. The lower row of Fig. 6 
essentially shows the situation for flat bands at the surface. 
These 1{E, k) images show the most ‘upshifted’ band struc¬ 
tures we have observed. These experimentally measured data 
differ by at most 20 meV from the values resulting from the 
time-infinite extrapolation shown as the right-most symbols 
for each material type in Fig. 5. 

Comparison of the upper and lower row of 1{E, k) images 
in Fig. 6 brings the differing alignments of the bands and 
the Fermi level sharply into focus. For the n-type compounds 
(Bi 2 Se 3 and Bi 2 Te 2 Se), the bands move such that the Fermi 
level goes from being well into the conduction bands (BB, 
see the top-left sketch of the electronic structure), to only 
just grazing them (flat-band, see the bottom-left sketch) IMl . 
Forthep-type compounds (BSTS2, BSTS1.5, BSTS1.46), the 
bands change such that the Dirac point moves from in excess 
of 250 meV binding energy (with the conduction bands clearly 
cutting the Fermi level, as shown in the top-right sketch of the 
band alignment), to Ejj lying above Ep, such that the valence 
band maximum is close to the Fermi level in the flat-band sit¬ 
uation (sketched in the lower-right schematic). 

Consequently, Fig. 6 is a nice illustration that once an 
effective flat-band-like situation can be reached, the Fermi 
level is either close to the bottom of the bulk CB or just above 
the top of the bulk valence band, as would be expected from 
straightforward semiconductor physics considerations. 

IV. CONCLUDING REMARKS 

All in all, our results offer new insight into an ongoing topic 
of great current interest, namely the optimization of the energy 
position of the topological surface states with respect to the 
Fermi level and with respect to the bulk bands in three dimen¬ 
sional topological insulators. Effective optimization strategies 
can be viewed as the first step towards future, application- 
oriented research. The new results, combined with the ad¬ 
vantage of hindsight, enable us to arrive at an interpretative 
framework into which both past experimental work should be 
placed and the context in which future studies can be viewed: 

(a) Fast downward band bending in UHV. From an 
ARPES point of view, the tunability of the energy position 
of the TSS cannot be solely accomplished by the bulk sto¬ 
ichiometry. 3D topological insulator compounds with bulk 


stoichiometries chosen so as to give real bulk insulating trans¬ 
port characteristics (e.g. BSTS1.5 and BSTS1.46) are subject 
to adsorbate induced band bending at their surfaces just as 
their binary counterparts lUMIllSllllEl. Adsorbate in¬ 
duced band bending in quaternary compounds evolves quickly 
after cleavage and details of the band structure measured us¬ 
ing ARPES depend as much if not more on the time interval 
between cleavage and measurement, than on the composition 
itself When the downward band bending saturates, a well- 
resolved conduction band is seen to cross the Fermi level at 
low temperature, thus although the bulk may be insulating, the 
surface region possesses more than only the TSS conduction 
channel. 

(b) Strong effect of illumination. The effect of a photon 
beam on the surface electronic structure of Bi 2 Se 3 , ternary 
and quaternary Bi-based TI compounds has been investigated 
and found to be considerable for fluxes typical at 3'^'* gen¬ 
eration synchrotron radiation sources, given an adsorbate- 
induced downward band bent situation as the starting point be¬ 
fore illumination. The direct and significant effect of the pho¬ 
ton beam underlines that great care needs to be taken when us¬ 
ing ARPES to infer quantitative trends in the energy position 
of the various features in the electronic band structure. Not 
only do the lifetime of the cleave and the background pres¬ 
sure (strictly also residual gas chemical composition) need to 
be logged, but also the photon flux and the irradiation time 
of the sample location under investigation. These findings 
may enable the generation of quantitative agreement between 
previous studies ll20l 1311 l34l ij?! . Interestingly, the value of 
photon flux used to effectively cancel the adsorbate-induced, 
downward band bending is comparable to that on offer in reg¬ 
ular laboratory laser and LED sources [1.3 x 10^^ photons/(s 
m2)]. 

(c) Combination of surface pbotovoltage and desorption 
phenomena. The marked effects of high-flux VUV illumina¬ 
tion on the energy landscape of band-bent Bi-based 3D TIs 
are attributed to the combined effects of surface photovoltage 
and photon-stimulated desorption phenomena. The combi¬ 
nation of an electronic process (SPV) and a mass-transport- 
related process (desorption) fits the dual ‘fast’ and ‘slow’ 
(twenty times slower than the ‘fast’) characteristic timescales 
observed for the photon-beam-related upshift of the electronic 
states. In Appendix 1, raising/re-lowering the sample tem¬ 
perature is shown to reversibly decrease/increase the down¬ 
ward band bending across the whole sample, and is due to 
temperature-induced desorption/re-adsorption at the surface. 
Appendix 2 shows that on a qualitative level, the upward en¬ 
ergy shift of key features of the electronic structure that ac¬ 
companies high-flux induced band flattening does not depend 
on the microscopic origin of the band bending, and also that 
the efficacy of this band-flattening is reduced at elevated tem¬ 
perature. Both these observations on the data presented in Ap¬ 
pendix 2 are in keeping with the importance of surface photo¬ 
voltage in the observed phenomena. 

(d) High-flux / high-energy illumination yields flat-band 
conditions at the surface. As the onset of downward band 
bending on adatom adsorption (even in UHV) is very rapid, 
using a standard ARPES experiment to reliably determine 
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the binding energies (relative to the Fermi level) of the Dirac 
point, as well as the valence and conduction bands in the near¬ 
surface region is not practicable. In contrast, extrapolation 
of the energies determined in high-flux, ARPES experiments 
conducted at higher photon energies is able to yield reason¬ 
able values for the relative alignment of the key features of the 
electronic structure, due to the effective flattening of the bands 
at and near the surface under such conditions. The higher is 
the flux X hi/ product, the closer is the extrapolated energy 
value to flat-band conditions. 

(e) Reliable ARPES determination of n- or p-type char¬ 
acter at low temperatures. Use of extrapolated high-flux 
(and/or high-h:^) binding energies (i.e. flat-band conditions) 
shows Bi 2 Se 3 and Bi 2 Te 2 Se to be n-type. On the other hand, 
BiSbTeSe 2 , and the two high-resistivity BSTS compositions 
studied [Bii, 5 Sbo, 5 Tei. 7 Sei .3 and Bii, 46 Sbo, 54 Tei. 7 Sei. 3 ] are 
p-type materials. 

Taken together, points (a)-(e) show how a full understand¬ 
ing of the energetics of the ARPES data from Bi-based 3D 
topological insulators can be arrived at, once explicit con¬ 
sideration is taken of the fact that illumination is an integral 
part of the experiment. Not only can high-flux ARPES data 
be used for a reliable determination of the character of the 
majority carriers in the TI system, but they also point to 
the possibilities of illumination-based, in-situ tuning of the 
key features of the electronic structure. Looking forwards, 
judicious use of super-band-gap illumination at a flux level 
readily achievable with laboratory light sources could be used 
as a valuable and reversible tool to decide whether additional 
topologically trivial surface-confined bulk states are at the 
chemical potential, contributing to transport, or whether 
future device-related experiments can be carried out using 
bulk crystals in which only the topological surface states 
count. 

METHODS 
Sample growth 

Crystals were grown in Amsterdam using the Bridg¬ 
man technique. High purity elements (Bi [99.999 %], Sb 
[99.9999 %], Te [99.9999 %] and Se [99.9995 %]) were 
melted in evacuated, sealed quartz tubes at 850°C and 
allowed to mix for 24 hours before cooling. The cooling 
rate was 3°C per hour. Samples were cleaved in UHV 
(P < 5 X 10“^°mbar) and at temperatures between 16-38K. 

Angle resolved photoelectron spectroscopy 

ARPES experiments were performed at four different ex¬ 
perimental setups; 

(i) At the SIS-HRPES end-station of the Swiss Light Source 
with a Scienta R4000 hemispherical electron analyzer. The 
minimum sample temperature was 16K and the pressure dur¬ 
ing measurements was 5 x 10“^^ mbar. 

(ii) At the UE-125 NIM beamline of BESSY II using the 
IDEEAA end-station 1651 with a VG Scienta ArTOE 10k elec¬ 
tron analyzer. The minimum sample temperature was 20K and 
the pressure during measurements was better than 5 x 10“^° 
mbar. 


(hi) At the UE112-PGM-2a-l2 end-station of BESSY II 
with a Scienta R8000 hemispherical electron analyzer. The 
minimum sample temperature was 35K and the pressure dur¬ 
ing measurements was 2 x 10“^° mbar. 

(iv) In Amsterdam using a Scienta VUV5000 helium pho¬ 
ton source and a Scienta 2002 hemispherical electron ana¬ 
lyzer. The minimum sample temperature was 19K and the 
pressure during measurement was 2 x 10“^° mbar (this is a 
pressure of ultra-pure He gas, the base pressure without oper¬ 
ating the light source was 3 x 10“^^ mbar). 

Data relevant to the competition of adsorbate- and photo- 
induced shifts were measured with hi/ = 23 eV and hi/ = 27 
eV, while the polarization of the photon beam was linear 
horizontal. Exposure to high-energy light was done with 
hi/ = 70 eV and hi/ = 130 eV. After exposure to high-energy 
photons the 1{E, k) images were recorded with hi/ = 23 eV 
or hi/ = 27 eV. The photon flux was deliberately decreased 
from Fhigh (SIS and UEl 12 sources only) to Flow by detuning 
the undulator gap, while keeping the nominal monochromator 
energy and exit slit - and thus photon footprint at the sample 
- fixed. Relative values of the photon flux were deduced by 
comparing with current values from a standard photodiode. 
The energy position of the Fermi level was determined within 
5 meV by measurements of an evaporated Au thin film that 
was held in electrical contact with the sample. 
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APPENDIX 1: the effect of sample temperature on BSTS1.5 

The data in this first part of the appendix show that chang¬ 
ing the sample temperature provides a fully reversible ap¬ 
proach by which the electronic properties at the surface can 
be changed on a global spatial level. Fig. Ala presents the 
energy position of the Dirac point and the CB minimum of a 
bulk-insulating BSTS1.5 crystal as a function of raising the 
sample temperature. The photon flux used for the measure¬ 
ment is low enough that conventional, downward BB is ob¬ 
served when the sample is continuously irradiated in UHV at 
fixed T. Fig. Alb shows the time-dependence of the relevant 
band energies for a different BSTS1.5 sample while it ther- 
malises after rapid heating of the cryostat from 20K to 150K. 
It is clear from both panels that the effect of increasing sam¬ 
ple temperature is to decrease the downward BB, as has been 
reported by others lISTlI^ . 

Fig. Ala shows that the shift of the bands to lower binding 
energy can be reversed again, simply by reducing the temper¬ 
ature once more. The open circles for T=20 K are the result 
of such a re-cool process. The small green arrows in Fig. 
7a show that after re-cooling, the binding energies of both 
the Dirac point and the CB continue to increase with time. 
As the temporal evolution of the downward shift (to higher 
binding energy) on cooling and thermalizing essentially tracks 
the time-scale seen for a freshly-cleaved sample exposed to 
residual vacuum at low temperature, this T-dependent shift is 
most likely to be due to the reversible adsorption/desorption 
of residual gas molecules at the TI surface ll20ll^ . 

Fig. Ala shows that for BSTS1.5, this means for temper¬ 
atures above ca. 200K, the CB in the near-surface region 
is above Ep and thus depopulated. The l(E, k) image for 
BSTS1.5 shown in Fig. Ale illustrates that the surface band 
bending for this material at room temperature is insufficient 
to lead to occupation of the CB, even after 15h of residual gas 
exposure. Nevertheless, for samples of macroscopic thickness 
the conductivity at higher temperatures will be dominated by 
impurity-related states M which do not give rise to coherent, 
sharp bands in fc-space. Fig. Ale shows that on re-cooling to 
40K, the ARPES signature changes significantly and only 10 
minutes at low temperature are enough to cause the downward 
BB to increase by a further 100 meV, which once again brings 
the CB minimum below Ep, as seen in the right-hand panel 
of Eig. Ale. Unlike the flattening of the band bending via 
the photo-induced effects described in the main body of the 
paper, the effects of changing the temperature are not local, 
but global in character. This is the reason that experimental 
broadening of the l(E, k) images like that seen in Eig. 3 is not 
observed in the data of Eig. Ale. 

The data of Eig. A1 - ported to the world of transport ex¬ 
periments - mean that even in micron thin BSTS flakes in 
which the bulk contribution can be suppressed during a trans¬ 
port experiment OOl . a second, topologically trivial conduc¬ 
tion channel will open up in the SCR at low temperatures. As 
a matter of fact, (1) adatom adsorption on a freshly cleaved 
sample held at low temperature and (2) cooling of a previously 
adatom-decorated sample will both lead to a fast downward 
energy shift with respect to the Eermi level, and to the occu¬ 
pation of the conduction band within the space charge region 


near the surface. 

The time-scale for this cooling-induced shift to higher 
binding energy (see Eig. Ala and the right two panels of Eig. 
Ale) is a lot faster than the time-scale of the return to down¬ 
ward BB seen after switching off the high-flux illumination 
that leads to band flattening (see the right-hand, grey-backed 
panel of Eig. 3a). After cessation of high-flux illumination 
at low temperatures, the photo-induced upward energy shift 
decays very slowly and only partially, and this points to 
an important role for SPV in the band flattening due to the 
relatively long recombination lifetime of the electrons and 
holes created by the photon beam. In contrast, on re-cooling, 
the bands shift back completely to higher binding energy, and 
do so at a rate resembling the fast, initial changes seen on the 
adsorption of residual gas atoms during the first few minutes 
of vacuum exposure at low T on a freshly-cleaved sample 
(see, for example, the data of Eig. 2c). These dynamics point 
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FIG. A1: The effect of temperature on the electronic band struc¬ 
ture of Bii.sSbo.sTei.rSei.s (BSTS1.5) (a) T-dependent energy po¬ 
sition of the Dirac point and the conduction hand (CB) minimum on 
a temperature ramp from 20 to 300K (filled circles). After 18 hours 
at 300K, the data shown hy the open green symbol were measured 
prior to re-cooling of the sample to 20K on which the binding energy 
of both the Dirac point and the CB minimum quickly recover to the 
previous, low-T values (uppermost open symbols at T=20 K). Both 
features continue shifting to higher binding energy as time passes (in 
this case a further 60 min), signalled by the small green arrows, (b) 
Time-dependent energy position of the Dirac point (left axis) and the 
valence band (VB) maximum (right axis) when the cryostat temper¬ 
ature is quickly increased to 150K but before full thermalization of 
the sample is attained. The error bars in the value of binding energy 
are estimated to be ±5 meV for the valence band and range between 
±10 and ±15 meV for the Dirac point, (c) Near-Tp electronic band 
structure at 300K measured 10 min after cleavage (left), 15 hours af¬ 
ter cleavage (middle) and upon subsequent cooling over 45 mins to 
40K, followed by a thermalization period (waiting time) of lOmins 
(right panel). 
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to the mechanism of the fast, spatially global modihcation 
of the electronic structure observed on changing the sample 
temperature as being desorption/re-adsorption of residual gas 
atoms on increasing/decreasing the sample temperature. This 
is in good agreement with results reported for both binary and 
ternary Bi-based TI materials 

APPENDIX 2: the effect of sample temperature on BSTS1.5 
decorated with metal adatoms 

Here we present the effect of high-flux photons on different 
kinds of equivalently band-bent BSTS1.5 surfaces held at dif¬ 
ferent temperatures. As starting points, we achieved a compa¬ 
rable amount of downward band bending at room temperature 
(RT) and at 38K by evaporating 0.1 ML of Ag at 35K (blue 
data points) and 0.8 ML of Nb at RT (red data points). Details 
of the development of the electronic structure of the interfaces 
between metals and BSTSL5 surfaces will be reported else¬ 
where. Here it is sufficient to note that - since the total value of 
photon-induced energy shift depends on the amount of initial 
band bending - a comparable degree of downward band bend¬ 
ing is important in order to perform a meaningful comparison 
between the situation at 35K and at room temperature. The 
central panel of Fig. A2 shows that high-flux photons have 
practically no effect on the electronic band structure of band- 
bent BSTSL5 held at RT (red data points). This is different 
to the situation where the same photon fluence affects a sim¬ 
ilarly band-bent BSTSL5 sample which is now kept at 35K 
(blue data points). In the latter case, the energy of the Dirac 
point shifts to lower binding energies by around 70 meV af¬ 
ter 20 minutes of photon exposure. There are two important 
conclusions that can be made from these data. Firstly, the 
qualitative effect of high-flux photons on BSTSL5 does not 
depend on the origin of the initial band bending. The same 
trends are seen, irrespective of whether the initial BB is a re¬ 
sult of the adsorption of residual gas atoms from UHV (Figs. 
1 -4 in the main paper) or due to deliberate surface decoration 
with metals (Fig. A2). Secondly, increasing the temperature 
to room temperature can effectively cancel the effect of high- 
flux photon exposure. This would be in keeping with strongly 
increased thermionic emission over the space charge region 
at high temperature resulting in an enhanced bulk-to-surface 
recombination current ll46ll66llF7l . The SPV mechanism pro¬ 
vides a natural fit to both of these observations, since it de¬ 
pends on the sample temperature and the amount of initial 
band bending but not on whether it arises from residual gas 
adatoms or from intentional metal adatom evaporation onto 
the surface. 
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FIG. A2: The effect of high-flux photons on Bii.sSbo.sTei.TSei.a 
(BSTS1.5) hand-hent due to metal adatoms and held at different 
temperatures. (a,b) The electronic band structure of BSTS1.5 after 
the deposition of 0.1 ML Ag at 35K. The spectra are acquired after 
(a) 1 min. and (b) 20 min. exposure to high-flux 27 eV photons 
at 35K, at which temperature the images were also measured, (c) 
The evolution of the Dirac point binding energy for the same band- 
bent BSTS1.5 after exposure to high-flux photons at 35K (blue data 
points) and at RT (red data points). The error bars in the value of 
binding energy are estimated to be ±5 meV. (d,e) The electronic hand 
structure of BSTS1.5 after the deposition of 0.8 ML Nh at RT. The 
spectra are acquired at room temperature after (c) 1 min. and (d) 23 
min. exposure to high-flux photons at RT. 
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